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Elevated-Temperature  Creep  Buckling  of 
Thermoplastic-Matrix  Fiber  Composites 
under  Biaxial  Loading 


Y.  Nakajo*  and  S.  S.  Wang** 

Department  of  Mechanical  and  Industrial  Engineering 
University  of  Illinois 
Urbana-Champaign ,  IL  61801 

ABSTRACT:  An  analytical  study  has  been  conducted  to  investigate  the  elevated- 
temperature  creep  buckling  behavior  of  thermoplastic-matrix  composite  laminates  sub¬ 
jected  to  multiaxial  loading.  Anisotropic  time-temperature  dependent  viscoelastic  consti¬ 
tutive  equations  of  the  composite  are  constructed  first,  using  the  experimental  results  and 
the  modified  Prony  series  expansion  solution.  Creep  buckling  loads  and  associated  creep 
failure  times  are  determined  by  a  time-dependent  bifurcation  buckling  analysis  of  the  ther- 
mopl  istic  composite  laminates  under  general  loading  in  an  elevated-temperature  environ¬ 
ment.  Detailed  solutions  for  high-temperature  creep  buckling  and  associated  failure  mode 
shapes  of  the  AS4/J1  thermoplastic  composite  laminates  with  various  fiber  orientations 
are  obtained.  The  creep  buckling  mode  shape  changes  and  the  accompanying  failure  loads 
are  studied  for  three  commonly  encountered  loading  modes:  (1)  pure  shear,  (2)  biaxial 
compression  with  any  arbitrary  stress  biaxiality  ratios,  and  (3)  combined  axial  compres¬ 
sion  and  shear  loading. 


1.  INTRODUCTION 

High-temperature  thermoplastic-matrix  composites  are  currently  con¬ 
sidered  to  be  used  in  high-performance  structures  such  as  supersonic  air¬ 
craft.  The  advanced  structural  applications  require  accurate  analyses  and  clear 
understanding  of  the  failure  behavior  of  structural  elements  at  elevated  tempera¬ 
tures.  Among  various  structural  failure  problems  presently  concerned,  elevated- 
temperature  stability  of  thermoplastic-matrix  composite  panels  under  mu'tiaxial 
loading  is  of  significant  interest  not  only  because  of  its  intrinsic  complexities  but 
also  because  of  its  technical  importance. 
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The  objectives  of  the  present  study  are  to  introduce  a  proper  mathematical 
method  to  account  for  high-temperature  constitutive  properties  of  the  composites 
and  to  investigate  fundamental  creep  buckling  behavior  of  thermoplastic-matrix 
composite  panels  under  general  multiaxial  loading  in  elevated-temperature  envi¬ 
ronments. 

The  subject  of  creep  buckling  of  structures  has  been  studied  by  many  research¬ 
ers,  and  two  common  approaches  have  been  taken.  One  is  based  on  nonlinear 
power  law  creep  stress-strain  constitutive  equation,  representative  of  constitutive 
laws  for  most  metal  creep.  Hoff  [10]  firstly  proposes  a  rational  criterion  for  this 
type  of  creep  failure  by  introducing  a  critical  time  (i.e.,  buckling  time)  at  which 
the  initial  imperfection  grows  to  infinite  amplitude.  A  brief  history  of  this  ap¬ 
proach  is  introduced  in  Reference  [11].  The  other  kind  of  creep  buckling  analysis 
is  the  one  used  for  structures  with  hereditary-integral-type  viscoelastic  material 
constitutive  equations.  The  definition  of  the  critical  time  used  in  this  approach 
with  small  deflection  theory  gives  an  infinite-long  buckling  time,  which  ap¬ 
parently  contradicts  to  experimental  results  [12].  To  alleviate  this  problem,  for 
example,  Hilton  defines  [7]  the  creep  buckling  time  of  a  column  as  the  one  at 
which  the  deflection  is  of  such  a  magnitude  that  the  maximum  moment  in  the  col¬ 
umn  will  reach  the  ultimate  value  the  material  can  carry.  Freudental  [8]  deter¬ 
mines  the  critical  time  of  a  column  by  considering  the  divergence  of  successive 
approximations  for  the  lateral  deflection  in  terms  of  a  power  series,  involving  the 
applied  load  and  time.  The  problem  of  creep  buckling  of  a  linear  viscoelastic 
material  has  been  investigated  widely  from  the  view  point  of  growth  of  an  initial 
imperfection  [2,14,15,16].  Williams  [4]  and  Powell  [5]  define  the  creep  buckling 
time  of  a  viscoelastic  column  based  on  a  bifurcation  buckling  criterion.  Recently, 
Wilson  and  Vinson  [21]  extend  the  viscoelastic  linear  buckling  analysis  to  lami¬ 
nated  composite  columns  using  the  correspondence  principle  to  obtain  the  criti¬ 
cal  time.  It  is  well-known  that  a  direct  application  of  the  bifurcation-type  analysis 
for  creep  buckling  generally  gives  extremely  modest  estimation  of  the  creep 
buckling  time  [10],  as  compared  with  Hoffs  approach.  However,  each  of  the 
above  analyses  should  give  certain  insights  into  the  creep  buckling  failure  be¬ 
havior  of  high-temperature  structures  with  different  configurations. 

In  this  research,  a  combined  experimental  and  theoretical  approach  is  taken  to 
analyze  the  creep  buckling  problem  of  a  thermoplastic-matrix  composite  panel  at 
elevated  temperatures.  The  experimental  results  obtained  from  accompanying 
studies  [17,20.22]  are  used  for  material  considerations.  A  recently  developed 
method  with  the  aid  of  an  optimization  theory  is  introduced  first  to  describe  the 
long-term  anisotropic  viscoelastic  behavior  of  the  composite.  Accuracy  and  con¬ 
vergence  of  the  solutions  are  studied  in  detail.  The  well-known  transform  method 
is  then  used  to  studv  the  creep  instability  of  the  thermoplastic-matrix  composite 
panels  with  several  different  fiber  orientations.  Three  loading  modes  are  con¬ 
sidered  in  the  study:  (1)  pure  shear,  (2)  biaxial  compression  of  arbitrary  biaxial 
stress  ratios,  and  (3)  combined  axial  compression  and  shear  loading.  The  results 
obtained  in  this  study  provide  important  information  to  our  understanding  of  the 
elevated-temperature  instability  of  the  thermoplastic-matrix  composite  panel 
under  general  loading  modes. 
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2.  METHODS  OF  APPROACH 

2.1  Anisotropic  Viscoelastic  Constitutive  Equations 
for  Unidirectional  Composites 

In  the  development  of  analytical  methods  to  study  elevated-temperature  creep 
buckling  behavior  of  a  thermoplastic-matrix  composite,  time-temperature- 
dependent  constitutive  equations  of  the  composite  must  be  determined  first.  In 
this  study,  the  anisotropic  viscoelastic  creep  compliance  data  of  a  unidirectional 
thermoplastic  fiber  composite  obtained  in  the  related  studies  [17,20,22]  are  used 
for  construction  of  material  constitutive  equations.  The  master  curves  of  creep 
compliances  of  the  composite  are  determined  by  applying  the  time-temperature 
equivalence  principle.  The  transverse  compliance  J22  and  ine  shear  compliance 
J66  of  the  composite  are  directly  measured  in  the  elevated-temperature  creep  ex¬ 
periments.  The  time-temperature-dependent  Poisson’s  ratio  pl2  is  determined  by 
the  use  of  the  rule  of  mixtures  [1],  based  on  the  properties  of  neat  matrix  resin 
and  the  fibers  at  different  temperatures  and  by  properly  shifting  the  data  to  con¬ 
struct  the  master  curve.  The  longitudinal  modulus  is  assumed  to  be  temperature- 
and  time-independent. 

The  composite  viscoelastic  constitutive  equations  are  expressed  in  the  follow¬ 
ing  forms  of  Equations  (2-1-1)  and  (2-1-2),  for  the  thermoplastic  matrix  com¬ 
posite. 


n 

J„(T,t)  =  a°'(T)  +  £  Ci‘>(D  exp  l^(T)t]  (2-1-1) 

k=\ 

n 

v<AT,t)  =  D!°'(T)  +  £  D<*'(D  exp  [7,(rv]  (2-1-2) 

k  =  I 

where  D\^\  nk  and  7*  are  coefficients  dependent  on  the  temperature  and 
fiber  orientations. 

To  express  the  experimental  data  of  the  anisotropic  composite  properties  by 
proper  analytical  expressions,  a  sequential  unconstrained  minimization  tech¬ 
nique  (SUMT)  is  used  in  this  study.  Details  of  the  SUMT  is  given  in  Appendix 
1.  Curve  fitting  for  a  longer  time  range  requires  the  introduction  of  several  new 
features  in  the  solution  procedure,  including  (1)  use  of  sufficient  number  of  terms, 
(2)  introduction  of  additional  constraint  conditions,  and  (3)  selection  of  prooer 
initial  values.  The  use  of  sufficient  terms  in  Equation  (2-1-1)  is  most  essential 
because  the  time  range  which  one  term  can  cover  is  limited  only  to  a  few  decades 
at  most,  and  the  regions  near  to  both  ends  of  the  time  range  of  interest  need  suffi¬ 
cient  terms.  An  accurate  solution  in  the  desired  time  range  requires  many  terms 
but,  on  the  contrary,  most  of  the  terms  are  used  to  fit  the  end  portions  of  the  curve 
(see  Appendix  1).  And  very  often  the  resolutions  are  lost  because  of  the  limit  of 
the  precision  in  the  calculation. 
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Besides  the  original  constraint  for  pk(t)  due  to  physical  feasibility,  additional 
constraint  conditions  are  introduced  not  only  to  accelerate  the  convergence  but 
also  to  prevent  the  solution  from  having  the  oscillating  feature,  which  sometimes 
introduces  a  small  error.  Besides  these  additional  conditions,  the  use  of  many 
terms  itself  increases  the  number  of  constraint  conditions 

With  increase  in  number  of  terms,  proper  selection  of  initial  values  becomes 
difficult.  The  SUMT  method  is  employed  to  solve  the  problem  efficiently  with 
above  three  requirements.  The  procedure  has  demonstrated  to  minimize  the  dis¬ 
crepancies  between  experimental  data  and  modified  Prony  series  solution  under 
these  constraint  conditions. 

2.2  Viscoelastic  Composite  Laminate  Theory 

Owing  to  the  elevated  temperature  environment  considered,  viscoelastic  com¬ 
posite  laminate  theory  needs  to  be  used  to  analyze  the  creep  buckling  behavior  of 
thermoplastic-matrix  composites  with  various  ply  orientations.  Specifically,  the 
classical  lamination  theory  is  modified  to  include  the  viscoelastic  composite  con¬ 
stitutive  equations.  The  use  of  the  well-known  correspondence  principle  [4]  facil¬ 
itates  the  analytical  procedure  of  the  composite  creep  buckling  analysis.  In  the 
present  study,  only  symmetric  angle-ply  composite  laminates  are  considered.  As 
a  result,  in-plane  and  out-of-plane  displacements  are  decoupled  in  the  viscoelas¬ 
tic  composite  plate  equations  [3] . 

Based  on  the  correspondence  principle,  the  elastic  composite  laminate  theory 
forms  the  basis  of  the  viscoelastic  composite  analysis.  Consequently  one  obtains 
the  corresponding  elastic  governing  equation  for  the  deflection  of  the  composite 
plate: 


_  d4w  d4w  d4w  d4w 

Dxx  —  +  4 Dxs  —  ..v  -f  2 (Dxy  +  2 Dss)  „  ,  +  4 Dv 


dx4 


dx3dy 


dx2dv2 


'ys 


dxdy3 


d4w  d2w  d2w  d2w 

+  Dyy  7TT  -I-  Nx  ~Z~T  +  7TT  ~  2 Nxv  „  7  =  0 


dy4 


dx2 


y  dy2 


dx  dy 


(2-2-1) 


where 


Daa  —  2 


n  eh„ 

L 

k=  1 


z2E{akg)dz  (a.  0  =  x,y,s) 


F  _  En  cos4  0 

‘■-'XX  1 

1  -  ^12^2, 
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\  1  -  ^,2^21 
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in  which  6  is  the  off  axis  angle,  and  the  numeric  subscripts  of  the  material  proper¬ 
ties  denote  those  in  material  coordinates  (Figure  2-2-1). 

The  solution  for  the  above  governing  differential  equation  along  with  proper 
B.C.’s  of  a  given  composite  panel  can  be  determined  by  various  methods.  The  vis¬ 
coelastic  solution  can  then  be  derived,  using  the  correspondence  principle  later. 

2.3  Viscoelastic  Creep  Buckling  Analysis  of 
Thermoplastic-Matrix  Composite  Laminates 

In  this  study,  Galerkin’s  method  [3]  is  used  to  obtain  the  corresponding  elastic 
buckling  solution.  The  elastic  solution  is  then  converted  into  viscoelastic  solu¬ 
tions  by  the  correspondence  principle  which  is  technically  the  inverse  Laplace 
transformation  of  the  elastic  solution  after  replacing  the  compliances  Ju  with  sJu 
and  the  Poisson’s  ratio  v u  with  svy.  By  replacing  the  associated  material  proper¬ 
ties,  the  elastic  solution  becomes  identical  to  the  viscoelastic  solution  in  the  La¬ 
place  domain,  or  5-domain.  Hence,  inverse  Laplace  transformation  of  the  elastic 
solution  with  the  replaced  material  properties  gives  the  viscoelastic  solution  in 
the  real  time  domain.  The  critical  time  of  structural  instability  for  a  ther¬ 
moplastic  composite  laminate  with  th?  viscoelastic  properties  is  subsequently 
calculated  from  the'  solution  obtained  by  above  referenced  technique  based  on  a 
bifurcation  buckling  analyses  of  a  composite  laminate. 

2. 3.  a  THE  GALERKIN’S  METHOD 

The  lateral  deflection  was  for  a  simply-supported  retangular  composite  plate  is 
assumed  as 


=  E  E 

m  =  I  n  =  I 


Cmrr  Sin 


m-KX  .  mry 
—  sin  — ; — 


(2-3-1) 


According  to  the  well-known  Galerkin’s  method,  the  following  equations  can  be 
obtained: 
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Figure  2-2-1.  Asymmetric  composite  laminates  with  [0J02/0J0f  fiber  orientations  under 
general  inplane  compression  and  shear  loading. 


where 


Sij  =  ( i/LX)2DXS  +  (  j/LyYDy, 

Fu  =  (i/Lx)4Dxx  +  2(i/Lx)Hj/Ly)2(Dxy  +  2  Dss)  +  ( j/Ly)4Dyy 


To  solve  for  Nx,  Ny,  and  Nxy,  each  equation  has  to  contain  only  one  unknown.  For 
the  pure  shear  case,  the  unknown  is  Nxy  and  Nx  =  Ny  =  0.  For  the  case  of  com¬ 
bined  axial  compression  and  shear  or  the  biaxial  compression,  either  one  of  the 
two  loads  must  be  the  unknown  and  the  other  must  have  a  known  value.  In  this 
study,  Nx  is  chosen  to  be  the  unknown  for  the  latter  two  cases,  and  the  constant 
values  of  other  loads,  i.e.,  Ny  for  the  biaxial  compression  and  Nxy  for  the  com¬ 
bined  axial  compression  and  shear,  are  changed  parametrically.  By  using  the 
Galerkin’s  method,  the  governing  differential  equation  is  discretized  into  simulta¬ 
neous  equations.  For  the  assumed  deflection  Equation  (2-3-1),  the  set  of  simulta¬ 
neous  equations  is  further  divided  into  two  sets  of  simultaneous  equations  (see 
Appendix  5).  The  solution  obtained  from  Equation  (2-3-2a)  has  contribution  of 
two  components  in  the  assumed  deflection  which  are  expressed  by  m  =  1,  n  = 

1  and  m  =  2,  n  =  2  in  Equation  (2-3-1)  shown  in  Figure  2-3-1.  And  the  solution 
from  Equation  (2-3-2b)  has  the  ones  which  correspond  to  m  =  1,  n  =  2  and 
m  =  2,  n  =  1.  No  other  mode  coupling  exists  except  for  those  combinations. 
The  lowest  load  should  be  chosen  among  the  four  roots  for  each  buckling  load 
obtained  from  Equations  (2-3-2a)  and  (2-3-2b). 


C||  tin  Cjj  t»t* 


.  Buckling  shape  corresponding  to  each  term  in  Equation  (2-3-1). 
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Figure  2-3-2.  Elastic  buckling  mode  shape:  (a)  pure  shear  loading;  (b)  combined  axial  compression  and  shear  (t,\  =  0.5t,\). 
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2.3.b  INVERSE  LAPLACE  TRANSFORMATION  FOR 
VISCOELASTIC  CREEP  BUCKLING 

Once  the  elastic  solution  is  obtained,  the  use  of  the  correspondence  principle 
leads  to  the  viscoelastic  solution.  For  example,  expanded  explicit  form  of  Equa¬ 
tion  (2-3-2b)  for  angle-ply  laminates  in  VE  domain  are: 
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where  Bx  to  Bxn  are  time-dependent  coefficients.  The  Bx  for  example  has  the 
form: 
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in  which 
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Solutions  for  general  angle-ply  laminates  including  the  degenerated  cases  of  a 
unidirectional  composite  by  letting  6  =  0  and  a  cross-ply  laminate  are  shown  in 
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the  Appendices  2,  3  and  4.  The  inverse  Laplace  transformation,  Equation 
(2-3-6),  can  be  calculated  by  the  use  of  convolution  integral.  For  example,  the 
1 2222  term  has  the  form 


>2222 


mn  Ct~p 

vi2{u)v\2ip  ~  u)du  1  v2x{\')v2x (r  ~  p  ~  v )dvdp 

0  j  0 


!  i 


*^22(£)  Jl2((j  — 


r 


-</ 


J eeiv)  J —  T  —  q  —  r])drjdqdr 

(2-3-7) 


The  mode  shape  change  is  also  calculated  by  making  use  of  the  correspondence 
principle.  The  ratios  cn/c‘22  and  c,2/c2,  can  be  conveniently  used  to  represent  the 
buckling  modes.  For  example,  Cu/c22  for  the  angle-ply  laminate  is  calculated 
from  the  simultaneous  equations  (see  Appendix  5)  from  which  Equation  (2-3-2a) 
is  obtained  with  a  condition  for  nontrivial  solution  for  cu  and  c22  as 
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(2-3-9) 


Applying  the  correspondence  principle,  the  above  equation  is  transformed  into 
the  real  time  domain  as: 
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4"  ^2  1  (  *0022  ‘\  1  22  ) 
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=  0 


(2-3-10) 


where  A,  to  A»  are  constants  (see  Appendix  2).  The  equations  to  be  solved  for 
other  modes  are  shown  in  Appendix  5.  The  ratio  of  the  amplitude  <',,/c22  can  be 
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calculated  for  the  applied  load  by  substituting  the  buckling  time  (for  example  100 
minutes  in  the  present  analysis)  to  determine  the  parameters  Itjkl. 

3.  RESULTS  AND  DISCUSSION 

To  illustrate  the  fundamental  behavior  of  elevated-temperature  creep  buckling 
of  thermoplastic-matrix  fiber  composites,  analytical  solutions  have  been  obtained 
for  a  carbon  fiber  (AS4)  reinforced  polyamide  (Jl)  matrix  composite  under 
various  in-plane,  biaxial  compression  and/or  shear  loading.  The  AS4/J1  ther¬ 
moplastic  matrix  composite  considered  here  has  a  glass-transition  temperature 
Tg  =  293 °F  (145°C).  Anisotropic  viscoelastic  constitutive  properties  are  deter¬ 
mined  by  the  use  of  the  elevated-temperature  experimental  data.  The  elastic 
modulus  along  fiber  direction,  Eu,  is  assumed  to  be  time-  and  temperature 
independent  with  the  value  of  1.75  x  107  psi.  Creep  buckling  behavior  of  the 
unidirectional  AS4/J1  composite  is  studied  first  to  provide  a  general  reference  for 
the  subsequent  analysis  of  more  complex,  multilayered  angle-ply  and  cross-ply 
thermoplastic  composite  laminates. 

For  both  cases  of  the  unidirectional  composite  and  multilayered  laminates  the 
panels  have  dimensions,  Lx  =  10  inch,  and  Ly  —  10  inch,  and  a  uniform  ply 
thickness  h  =  0.2  inch  (Figure  2-2  1).  Three  multiaxial  loading  modes  are  con¬ 
sidered:  (a)  pure  shear  loading,  (b)  combined  axial  compression  and  shear,  and 
(c)  biaxial  compression  with  arbitrary  stress  biaxiality  ratios.  For  each  case,  a 
schematic  is  shown  in  the  respective  figure.  The  nominal  shear  stress  and  the 
nominal  compressive  stress  in  the  y  direction  are  changed  parametrically  in  the 
analyses  for  the  cases  (b)  and  (c).  respectively.  In  all  of  the  cases  studied,  creep 
buckling  loads,  critical  instability  times,  and  associated  buckling  mode  shapes  of 
the  composite  laminates  under  various  combinations  of  in-plane  biaxial  compres¬ 
sion  and  shear  are  determined  for  several  temperatures.  The  degenerated  solu¬ 
tion,  i.e.,  the  thermoelastic  buckling  solution,  for  each  case  is  also  determined 
and  included  here  for  the  reference  purpose. 

3.1  Analytical  Expression  of  Anisotropic 
Viscoelastic  Composite  Constitutive  Equations 

In  the  use  of  the  correspondence  principle  for  studying  the  thermoplastic  com¬ 
posite  constitutive  equations,  only  the  time  dependency  is  taken  into  account. 
The  temperature  effect  can  be  introduced  into  the  analyses  by  the  time-tempera¬ 
ture  equivalence  concept  [19 J.  Owing  to  the  fact  that  the  time  dependency  of  the 
composite  results  mainly  from  the  matrix  [22 J,  the  time-temperature  equivalence 
is  isotropic.  Therefore  the  shift  factor  of  the  matrix  can  be  used  to  relate  any  of 
the  solutions  obtained  in  the  present  analysis  at  145 °C  to  that  at  a  desired  temper¬ 
ature.  While  any  solution  at  a  given  temperature  can  be  shifted  to  that  at  a 
different  temperature,  the  expressions  of  the  viscoelastic  material  properties  in 
the  present  analysis  is  valid  only  for  the  time  range  of  five  decades. 

The  numbers  of  terms  used  in  the  modified  Prony  series  expansion  of  J22.  J,,„ 
and  v 1 2  are  determined  for  the  solutions  within  an  admissible  range  of  approx¬ 
imation.  The  present  method  can  theoretically  give  precise  fitting  for  a  long 
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range  by  using  increased  number  of  terms.  However,  owing  to  the  resolution  of 
numerical  results,  the  use  of  many  terms  (e.g.,  n  >  15)  may  now  work  out  for 
too  long  a  time  range.  At  present,  the  upper  limit  of  the  time  range  is  5  decades 
in  a  a  log  scale,  which  can  be  shifted  to  another  portion  of  the  master  curve.  But 
it  is  very  difficult  to  expand  to  the  range  of  more  than  5  decades.  Figure  3-1-5 
shows  the  difficulties  to  extend  the  ends  by  adding  additional  terms.  Sometimes 
close  values  of  the  design  variables  /3/s  lead  to  the  same  local  optimum  value  (see 
Table  3-1-1).  The  more  this  happens,  the  greater  the  advantage  of  starting  the  op¬ 
timization  using  many  terms  is  lost.  In  the  present  analyses  initial  values  are  de¬ 
termined  by  the  next  equation. 

For  /3/s, 


(3, =  1  -  g‘  -  h  (i  =  1 ,/!)  (3-1-1) 

where  0  <  g  <  1,  h  <  1.  For  C,s  (after  normalization), 

C()  =  1.2 

C,  =  -0.01  (/  =  1,/t)  (3-1-2) 

The  order  in  Tables  3-1-1  to  3-1-3  corresponds  to  the  order  of  the  initial  values  in 
Equation  (3-1-1). 

Convergence  of  the  series  solution  with  different  numbers  of  terms  can  be  seen 
from  Figure  3-1-5  which  shows  the  feature  of  each  normalized  term  for  represen¬ 
tative  (3,  values  in  Equation  (2-1-1).  To  extend  the  range  to  the  direction  of  t  =  0. 


Table  3-1-1.  Various  terms  in  modified  Prony  series  for  composite  creep 
compliance  J22af  T  =  293°F  (145°C). 


k 

C&Mpsi’) 

1 

-  1 .452100E  -  07 

-2.514420E  +  00 

2 

-  1 .146764E  -  07 

-  1.27771  IE +  00 

3 

-  9.236675E  -  08 

-  7.365715E  -  01 

4 

-  6.73441 7E-  08 

-  3.161 159E  -  01 

5 

-  5.990098E  -  08 

-  1  569684E  -  01 

6 

-  6.721 540E  -08 

-  8.364243E  -  02 

7 

-  8.490632E  -  08 

-  4.535500E  -  02 

8 

-  1 .027089E  -  07 

-  2.200023E  -  02 

9 

-  1.1 34783E- 07 

-  7.457791  E-  03 

10 

-  1  .  189835E  -  07 

-  6.627651  E- 04 

11 

-  1.202688E-07 

-  1  065326E  -  03 

12 

-  1 .221 522E  -  07 

-  1.602799E-04 

13 

-  1  227818E  -  07 

-  2.8691 80E  -  05 

14 

-  1  227750E  -  07 

-  2.8691 80E- 05 

15 

-  1.228800E-07 

-  2.8691 80E- 05 

Table  3-1-2.  Various  terms  in  modified  Prony  series  for  composite  creep 
compliance  J66  at  T  =  293°F  (145°C). 


k 

C£j  (Psr1) 

Mk 

1 

-2.00391  IE -07 

-  2.501 196E  + 00 

2 

-  1.658524E-07 

-  1.275838E  +  00 

3 

-  1 .4541 92E  -  07 

-  7.355021  E- 01 

4 

-  1. 342051 E -07 

-  3.151077E  -  01 

5 

-  1. 507481 E- 07 

-  1.553528E-01 

6 

-  1 .9001 40E  -  07 

-  8.061813E  -  02 

7 

-  2.4331 90E -07 

-  3.992683E  -  02 

8 

-  2.891804E  -  07 

-  1 .451749E  -  02 

9 

-  3.089745E  -  07 

-  3.267287E  -  03 

10 

-  3.123734E  -  07 

-  1. 37041 8E -03 

11 

-3.129629E  -07 

-  5.621973E  -  05 

12 

-  3.1 28881 E  -  07 

-  5.621 973E-  05 

13 

-  3.126838E  -07 

-  5.621 973E  -  05 

14 

-  3.124806E  -  07 

-  7.796962E  -  04 

15 

-  3.124339E  -  07 

-  1  155590E  -  03 

Table  3-1-3.  Various  terms  in  modified  Prony  series  for  Poisson’s  ratio  v12 

atT  =  293°F  (145°C). 


k 

Q(K) 

7* 

1 

- 1 .015074E  -02 

-  2.497946E  +  00 

2 

-  6.75561 5E- 03 

-  1.27491 6E  + 00 

3 

-  5.39631 3E- 03 

-  7.349936E  -  01 

4 

-3.69051  IE -03 

-  3.148159E  -  01 

5 

-  2.069948E  -  03 

-  1  554924E  -  01 

6 

-  1 .485730E  -  03 

-  8.2831 69E- 02 

7 

-4.751751E  -03 

-  4.716623E  -  02 

8 

-  8.328338E  -03 

-  2. 804054 E  -  02 

9 

-  1 .099923E  -  02 

-  1.58641  IE  -02 

10 

-  1. 283541 E -02 

-  5.166782E  -  03 

11 

-  1 .510358E  -  02 

-  4.046428E  -  05 

12 

-  1 .597748E  -  02 

-  3.410993E  -  04 

13 

- 1 .621792E  -02 

-  4.046428E  -  05 
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T  ‘  T0 

Fit  ure  3-1-1.  Shift  factor  vs  temperature  for  AS4/J1  thermoplastic-matrix  composite  [ refer- 
en  ■ e  temperature  =  293°F  (145°C)]. 


the  use  of  many  terms  with  very  small  (3,  values  is  necessary.  To  extend  the  range 
to  the  direction  of  t  =  oo,  additional  terms  with  very  close  (3,  values  are  neces¬ 
sary.  In  both  cases,  the  resolution  is  easily  lost  by  numerical  errors  during 
ca  culations.  Also  in  this  figure,  it  is  apparent  that  a  few  terms  used  can  cover 
only  2  or  3  decades  over  a  log  time  scale  at  most,  which  is  independent  of  the 
scheme  of  optimization  or  accuracy  of  calculation. 

Figure  3-1-1  shows  the  relation  between  the  temperature  and  the  shift  factor  by 


,  CAT  -  To) 

°g  °T  ~  C2  +  T  -  T0 


(3-1-3) 


which  is  sometimes  referred  to  as  the  WFL  equation  [19].  Values  of  the  C,  and 
C2  ire  experimentally  determined  constants,  Fis  a  temperature  and  T0  is  a  refer¬ 
ence  temperature.  Conventionally  Cx  and  C2  are  obtainable  by  a  graphical 
method  as  a  gradient  and  an  intercept  in  the  ( T  —  F0)/log  a,  vs.  T  —  T0  curve. 
In  the  present  data  reduction,  a  least-square  method  is  employed  to  obtain  C,  and 
C2.  When  a  glass  transition  temperature  F*  is  used  for  T0 ,  it  is  known  that  C,  and 
C2  give  approximately  consistent  values  regardless  of  the  material  for  the  temper¬ 
ature  range  T  =  F„  —  50 °C  to  T  =  TH  4-  50 °C.  The  C,  and  C2  values  obtained 
from  our  experimental  results  differ  from  those  consistent  values  [19]. 


J?*)  fo°i 


oo 


Figure  3-1-3.  High-temperature  creep  experimental  data  and  the  modified  Prony  series 
s  'lution  for  AS4/J1  thermoplastic  composite  at  T  =  293°F  (145°C). 
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.2. (•{*«•  (.«•£*•■  2. *•£«••  4.HE<M  S . fit «<• 

log  c 

Figure  3-1-5.  Each  normalized  term  for  representative  (1,  values  in  Equation  (2-1-1). 


3.2  Creep  Buckling  of  Unidirectional  Thermoplastic-Matrix  Composites 

j  La  PURE  SHEAR  LOADING 

As  a  reference  for  studying  the  elevated-temperature  creep  buckling  of  more 
complex  composite  laminates,  thermal  elastic  buckling  of  a  unidirectional  com¬ 
posite  panel,  i.e.,  rf2  =  t?2  (0.  Tg )  is  introduced  first.  The  rf2  in  Figure  3-2-1 
is  the  elastic  bifurcation  buckling  load  calculated  with  the  composite  ther¬ 
moelastic  properties  at  the  glass  transition  temperature  T  =  145 °C.  The  value  of 
rf2  is  19.36  ksi.  The  relationship  between  creep  stress  and  creep  buckling  failure 
time  of  the  composite  at  various  temperatures  is  shown  in  the  figure.  The  shift 
factor  aT  is  used  for  calculating  buckling  failure  at  various  T s.  The  results  at  the 
glass  transition  temperature  Tg  ( =  293  °F)  are  also  presented  for  comparison 
purposes.  Above  the  glass  transition  temperature,  the  creep  buckling  failure 
becomes  apparently  more  sensitive  to  the  applied  load.  In  other  words,  a  small 
change  in  the  applied  load  leads  to  a  large  change  in  the  creep  buckling  time  at 
T  >  Tg.  In  the  case  of  an  elastic  bifurcation  buckling,  the  critical  load  can  be 
un  quely  determined  for  a  given  composite  structure.  However,  in  the  present 
cr  ep  buckling  studies,  creep  buckling  loads  depend  upon  the  time-dependent 
m  iterial  properties.  For  more  complex  loading  modes,  the  creep  buckling  results 
arc  shown  in  Figures  3-2-1  to  3-2-7  for  both  unidirectional  and  multilayered 
angle-ply  and  cross-ply  thermoplastic  composite  laminates. 

3.2.b  BIAXIAL  COMPRESSION 

In  Figure  3-2-2  the  relationship  between  the  biaxial  creep  stresses  aT  (t.T) 
and  a f  (t,T)  and  the  associated  creep  failure  times  t  is  shown.  [The  buckling 


*12  (C,T)/r12 


I  .00E  *00 


* 


9.50E-01 


9.00E-01 


0.30E-91 


l.E-91  1  .  E  *00  I.E+B1  1  .E  *92  1  .E«93 


t/t* 

Figure  3-2-1.  Creep  buckling  failure  of  AS4/J1  thermoplastic-matrix  composite  panel  under 
pure  shear  loading  at  various  temperatures  (t '  =  1  min.). 
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Figure  3-2-4.  Creep  buckling  failure  of  AS4/J1  thermoplastic-matrix  composite  angle  ply 
under  pure  shear  loading  at  various  temperatures  (t’  -  7  min.). 
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Figure  3-2-6.  Creep  buckling  failure  of  AS4/J1  thermoplastic-matrix  composite  angle  ply 
under  combined  compression  and  shear  at  T  =  29 3°F  (14 5°C). 
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load  <r,“  ( t,T )  is  normalized  with  respect  to  of  which  is  the  thermal  elastic 
buckling  toad  (<rf  ==  6.974  ksi)  at  T  =  145 °C  under  uniaxial  compression.  The 
in  Figure  3-2-2  is  the  corresponding  reference  load  of  the  same  composite 
panel  under  uniaxial  compression  along  the  x2  axis  at  T  =  145  °C  and  t  =  100 
min.]  The  top  solid  line  in  Figure  3-2-2  is  the  creep  buckling  solution  for  o f  = 
0,  i.e.,  a  uniaxial  compressive  buckling  in  xx  direction.  The  essential  features  of 
the  biaxial  creep  buckling  in  Figure  3-2-2  are  similar  to  those  in  a  pure  shear 
buckling  case.  At  a  given  time  t  the  introduction  of  oT  reduces  the  critical  load 
a ?  ( t,T )  and  changes  the  mode  shape.  The  cusps  in  Figure  3-2-2  made  by  two 
solutions  corresponding  to  different  mode  shapes  have  significant  physical  impor¬ 
tance.  One  is  that  they  are  the  points  at  which  both  of  the  two  mode  shapes  are 
admissible.  The  other  is  that  sensitivity  of  the  buckling  time  to  the  applied  load 
a T  ( t,T )  changes  near  these  points.  The  most  significant  difference  from  the 
elastic  buckling  is  that  in  the  case  of  elastic  buckling,  the  mode  shape  change  cor¬ 
responds  to  the  change  in  specific  buckling  characteristics  of  the  structure.  But  in 
the  case  of  creep  buckling,  the  mode  shape  change  occurs  along  with  the  creep 
buckling  time. 

3.2.  c  COMBINED  AXIAL  COMPRESSION  AND  SHEAR  LOADING 

In  Figure  3-2-3,  the  relationship  between  the  combined  creep  stresses  oT 
( t,T )  and  t?2  ( tj )  and  the  associated  creep  failure  time  t  is  shown.  (Thermal 
elastic  buckling  loads  af  are  defined  as  that  in  the  case  3.2. b.)  The  top  solid  line 
is  identical  to  that  in  Figure  3-2-2,  i.e.,  the  uniaxial  compressive  creep  buckling 
along  the*i  axis.  Again,  introduction  of  additional  stress  rf2  reduces  the  critical 
applied  load  at  any  buckling  time.  The  shear  stress  tT2  is  the  reference  shear 
stress  which  alone  produces  a  100  min.  creep  buckling  time.  In  this  loading  case, 
the  reduction  is  not  uniform  even  without  a  mode  shape  change.  A  parametric 
study  of  the  combined  loading  shows  that  the  solution  is  much  affected  by  the 
presence  of  the  shear  stress  than  that  of  transverse  compression. 

3.3  Creep  Buckling  of  Thermoplastic  Composite  Laminates 

3. 3.  a  ANGLE-PLY  COMPOSITE  LAMINATES 

i)  Pure  Shear  Loading 

The  unique  feature  of  the  solution  for  an  angle-ply  thermoplastic  AS4/J1  com¬ 
posite  laminate  is  similar  to  that  of  the  unidirectional  composite  panel.  (/*,  in 
Figure  3-2-4  is  23.97  ksi  for  the  angle-ply  composite  laminate  T  =  145 °C.)  A 
comparison  with  Figure  3-2-1  shows  that  the  creep  buckling  critical  time  is  sig¬ 
nificantly  extended  for  a  given  applied  load.  It  is  hard  to  observe  directly  from 
these  figures  because  they  are  normalized  by  different  factors.  However,  the  creep 
buckling  property  or  sensitivity  of  the  creep  buckling  time  to  the  applied  load  is 
also  affected  by  the  composite  lamination  variables. 

ii)  Biaxial  Compression 

The  same  changes  can  be  seen  in  Figure  3-2-5.  (a*  in  Figure  3-2-5  is  11.20 
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ksi.)  The  mode  shape  change  observed  in  the  case  of  the  AS4/J1  unidirectional 
composite  does  not  appear  for  the  angle-ply  composite  laminate.  Figure  3-2-8 
also  shows  that  mode  shape  does  not  change  in  a  [45/  —  45/  —  45/45]  composite 
laminate. 

iii)  Combined  Axial  Compression  and  Shear 
Creep  budding  resistance  is  improved  similarly  as  the  above  two  cases.  The  re¬ 
duction  of  the  critical  load  due  to  presence  of  additional  shear  stress  at  a  given 
time  is  much  more  uniform  compared  to  that  of  the  unidirectional  composite. 
Thus,  the  creep  buckling  time  of  the  angle-ply  laminate  will  not  reduce  so  drasti¬ 
cally  as  the  case  of  the  unidirectional  composite  by  the  presence  of  the  shear 
stress. 

3. 3.  b.  CROSS-PLY  COMPOSITE  LAMINATES 
In  the  case  of  cross-ply  thermoplastic  laminates,  solutions  for  pure  shear  load¬ 
ing  and  combined  axial  compression  and  shear  loading  are  exactly  the  same  as 
those  of  the  unidirectional  composite.  The  reason  is  that  the  flexural  rigidities  in 
Equation  (2-3-2a)  are  not  changed  in  the  case  of  Lx  =  Ly  but  only  those  in  Equa¬ 
tion  (2-3-2b)  are  affected.  Consequently  only  the  solutions  which  have  the  ad¬ 
missible  roots  form  the  Equation  (2-3-2b)  have  different  values.  It  should  be 
noted  that  in  the  case  of  general  rectangular  laminates  which  have  different  Lx  and 
Ly ,  all  of  the  three  solutions  are  different  from  those  of  unidirectional  com¬ 
posites. 


0.00t«00  4 .  30C  +01  9 . 00E  *01  1.3SC«42  !.09f*0? 

9 

Figure  3-2-8.  Applied  load  which  causes  the  100  min.  buckling  time  vs.  angle  6  of  the  com¬ 
posite  angle  ply  [±0)s  under  biaxial  compression  at  T  =  293°F  (145°C). 
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Figure  3-2-7  shows  that  the  admissible  solution  from  Equation  (2-3-2b)  repre¬ 
sents  much  improved  creep  budding  resistance  so  that  the  failure  life  increases 
significantly.  Also,  transition  of  the  mode  shape  is  much  suppressed  as  compared 
with  those  of  the  unidirectional  composites. 

3.4  Mode  Shape 

Figure  3-2-8  shows  variations  of  the  applied  load  which  causes  the  100  min. 
buckling  time  in  the  biaxial  compressive  loading.  The  dotted  line  shows  the  ad¬ 
missible  border  between  the  two  possible  modes.  Figure  3-2-9  shows  the  similar 
curves  with  Figure  3-2-8  in  the  case  of  the  combined  shear  and  compressive 
loading.  The  dotted  lines  show  the  borders  along  which  the  admissible  solution 
changes  between  one  from  Equation  (2-3-2a)  and  another  from  Equation  (2-3- 
2b)  that  means  the  change  of  mode  shape.  The  off-axis  angle  0  =  0°  and  45° 
correspond  to  a  unidirectional  composite  and  angle-ply  laminate  respectively 
which  were  discussed  previously. 

Figure  3-2-10  and  3-2-11  show  how  the  large  fraction  of  the  shear  component 
can  change  the  mode  shapes  drastically  in  the  case  of  the  combined  axial  com¬ 
pression  and  shear  loading  that  is  not  seen  in  the  other  loading  case. 

4.  CONCLUSIONS 

Based  on  the  analytical  methods  developed  for  studying  the  creep  buckling 
failure  of  fiber  composite  laminates  and  the  results  obtained  for  the  AS4/J1 


# 


Figure  3-2-9.  Applied  load  which  causes  the  100  min.  buckling  time  vs.  angle  6  of  the  com¬ 
posite  angle  ply  [±6]s  under  combined  axial  compression  and  shear  at  T  =  293°F 
(145°C). 


Cll/P22 


Figure  3-2-11.  Buckling  mode  of  AS4/J1  thermoplastic-matrix  composite  angle  ply  under 
combined  axial  compression  and  shear  at  T  =  293°F  (145°C). 
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thermoplastic-matrix  composite  panels  under  various  loading  modes  in  the 

elevated-temperature  environment,  the  following  conclusions  can  be  made  in  this 

study: 

1.  The  elevated-temperature  anisotropic  viscoelastic  constitutive  equations  of 
the  thermoplastic-matrix  composite  can  be  constructed,  using  the  modified 
Prony  series  expansion  with  the  aid  of  a  sequential  unconstrained  minimiza¬ 
tion  technique.  The  series  solution  for  the  composite  time-temperature- 
dependent  properties  matches  well  with  the  experimental  data  in  a  range  of 
five  decade  along  the  logarithmic  scale. 

2.  The  elevated-temperature  viscoelastic  creep  buckling  of  the  composite  panels 
is  analyzed  successfully  by  the  use  of  the  well-known  Galerkin’s  method  and 
the  inverse  Laplace  transformation  technique.  At  a  specific  time  creep  buck¬ 
ling  loads  of  a  thermoplastic  composite  laminate  panel  are  determined  at 
various  temperatures  and  under  various  combinations  of  biaxial  loading 
modes. 

3.  The  creep  buckling  failure  time  of  a  thermoplastic  composite  panel  under  a 
given  load  (e.g.,  uniaxial,  shear  or  biaxial)  is  determined  by  a  time-dependent 
bifurcation  analysis  of  the  characteristic  equations  in  conjuction  with  the  in¬ 
verse  Laplace  transformation. 

4.  At  a  given  temperature,  the  shear  creep  buckling  load  of  the  unidirectional 
AS4/J1  thermoplastic-matrix  composite  appears  to  decrease  more  rapidly 
with  time  than  that  of  an  angle-ply  [±45°],  composite  laminate.  As  antici¬ 
pated,  the  symmetric  angle-ply  [±45°]s  laminate  has  a  better  resistance  to 
elevated  temperature  shear  creep  deformation  and  failure. 

5.  The  elevated-temperature  creep  buckling  resistance  of  a  composite  laminate 
panel  is  strongly  affected  by  the  biaxial  compressive  loading.  The  introduc¬ 
tion  of  of  in  a  compressively  loaded  (aT)  composite  plate  significantly 
reduces  the  axial  creep  buckling  resistance  in  all  three  cases  studied  (i.e., 
pure  shear,  biaxial  compression,  and  combined  shear  and  compression).  It 
also  changes  the  buckling  shape  modes  in  the  cases  of  uniaxial  composite  and 
cross-ply  composite  laminates.  However,  this  phenomenon  is  not  observed  in 
the  [±45°],  laminate  at  any  temperatures. 

6.  The  elevated-temperature  creep  buckling  strength  of  a  composite  laminate 
panel  is  also  affected  by  the  combined  axial  compression  and  shear  loading. 
The  additional  shear  introduced  to  the  axially  loaded  composite  significantly 
reduces  the  composite  panel  load  bearing  capacity.  However,  in  both  the  uni¬ 
directional  composite  and  the  [  ±45°]s  composite,  mode  shape  changes  at  any 
temperatures  and  load  combinations  are  not  observed. 

7.  In  biaxially  compressed  thermoplastic  angle-ply  composite  laminates  (i.e., 
of  0 ,  of  ^0),  fiber  orientations  affect  not  only  the  creep  buckling  load 
level  but  also  the  mode  shape  associated  with  the  failure.  For  example,  in  a 
[6/  —  6/  —  6/6 1  AS4/J1  composite  laminate  at  145°C  subjected  various  combi¬ 
nations  of  of/of,  the  dominant  mode  shape  during  creep  buckling  of  the 
laminate  with  6  <  33°  ism  =  1  and  n  =  1,  whereas  the  6  >  33°  the  buck¬ 
ling  mode  shape  is  m  —  1,  n  =  2  under  any  combination  of  cr *  and  of . 
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8.  In  the  case  of  combined  compression  and  shear,  fiber  orientation  does  not  af¬ 
fect  the  mode  shape  change  significantly  but  the  load  level.  Due  to  the  pres¬ 
ence  of  the  shear,  symmetry  of  the  characteristics  in  which  the  load  level  is 
affected  by  the  fiber  orientation  along  6  =  90°  is  lost. 

APPENDIX  1:  SEQUENTIAL  UNCONSTRAINED 
MINIMIZATION  TECHNIQUE  (SUMT) 

To  express  the  experimental  data  of  the  anisotropic  composite  properties  by 
proper  analytical  expressions,  a  sequential  unconstrained  minimization  technique 
[9]  which  is  one  of  the  commonly  used  optimization  techniques  is  employed  in¬ 
stead  of  the  conventional  Prony  expansion  method  [18].  In  the  present  study,  15 
terms  are  used  for  the  creep  compliances  and  13  terms  for  the  time-dependent 
Poisson’s  ratio.  The  number  of  the  terms  used  is  determined  by  two  factors  which 
conflict  with  each  other.  An  accurate  solution  in  the  desired  time  range  requires 
many  terms  but,  on  the  contrary,  most  of  the  terms  are  used  to  fit  the  end  portions 
of  the  curve.  And  very  often  the  resolutions  are  lost  because  of  the  limit  of  the 
precision  in  the  calculation. 

The  original  objective  function  /  for  the  present  problem  is 

N 

/  =  E  I /%,«,)  -  Fmp(t<)i‘  (A-l-1) 

(=  I 


where 

Fex(ti)  =  experimental  data  for  the  compliances  or  the  Poisson’s  ratio  at  /  = 

ti 

Fmp(t,)  =  the  compliance  function  or  Poisson’s  ratio  calculated  by  the  modified 
Prony  series  solution 

The  constraint  conditions  should  be  expressed  as  follows  to  be  used  with  SUMT 
method: 


gi  >0  (/  =  1  ,NC  ;  Nc:  Number  of  constraint  conditions) 

(A-l-2) 


For  the  compliances  for  example. 


0k(T)  >  0  (A-l-3a) 

where  /3*  is  the  new  variable  introduced  for  convenience  instead  of  nk.  The  rela¬ 
tions  between  them  are 


(3k(T)  =  exp  [^(T)] 


(A-l-4) 
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Besides  this  essential  constraint  condition,  two  additional  constraint  conditions 
are  necessary  due  to  the  reasons  mentioned  previously: 


>o  (k  =  \,n) 


(A-l-3b) 


1  -  0k(T )  >0  (k  =  l,#i) 


(A-l-3c) 


Similar  conditions  are  necessary  for  the  constants  in  the  expression  for  Poisson’s 
ratio. 

A  number  of  these  constraint  conditions  can  be  automatically  satisfied  by  mod¬ 
ifying  the  objective  function  as 


u  =  f  +  Yt  (r</£<) 


(A-l-5) 


The  summation  term  is  the  penalty  term  introduced  instead  of  the  constraint  con¬ 
ditions.  By  reducing  r,  gradually  to  zero  during  the  iteration  of  the  conventional 
steepest  descent  method  [9],  the  optimum  solution  for  u  car  be  obtained  which 
is  actually  the  same  solution  for  /with  all  of  the  constraint  conditions. 

The  SUMT  method  resolves  the  complicacy  to  handle  a  number  of  constraint 
conditions  but  still  difficulties  remain  in  selecting  proper  initial  values  and  in  re¬ 
ducing  the  magnitude  of  the  coefficients  of  the  penalty  terms  appropriately  which 
are  at  present  solved  by  trial  and  error. 

APPENDIX  2 

The  explicit  forms  for  the  flexural  rigidities  for  the  angle-ply  laminates  are: 

ry  _ _ _ [E22  sin4  0  +  2\v2iEu  +  2(1  —  ri2r2,)G66} 

IX  3(1  -  pi2v2l)  sin2  0  cos2  0  +  Eu  cos4  0J 

n  _  2r3  [£,,  sin4  0  +  2jr2I£1I  +  2(1  -  r12r21)G66} 

yy  3(1  -  vl2v2i )  sin2  6  cos2  6  +  E22  cos4  0J 

n  _  2* 3  [v2iEn(oLU4e  +  cos4  0)  +  t v2lEu  -  4 

3(1  -  vl2v2l)  (1  -  r12r2, )G66}  sin2  0  cos2  0J 


= 


[(£22  +  Eu  —  2r2j)  sin2  0  cos2  0  -f  G6 


3(1  -  vi2v2l)  (cos2  0  -  sin2  0)2] 


t3  [E22  cos  0  sin3  0  —  £n  cos3  0  sin  0  + 

- +  2(1  -  vx2v2i  )G66 }  sin*  cos8 

VxlVlx)  (sin2  0  -  cos2  0)] 

f 3  [E22  sin  0  cos3  0  —  Eu  sin3  0  cos  0  -l- 

- -  (^21  Eu  +  2(1  -  r12r2,  )G66 1  sin9  cos9 

Vl2Pu)  (sin2  0  -  cos2  0)] 
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APPENDIX  3 

The  expanded  form  of  Equation  (2-3-2a)  for  the  angle-ply  laminates  to  solve 
for  the  applied  loads  is: 

(4)022  —  2/u22  +  ^2222)^  xy  +  (20/3)7T2  +  2C3)Mi(/00,2  —  Am)  ~  A2l(  ^0022 
—  A122)  "1"  A3(Iq122  —  A222)  +  /44(/0021  —  2/  1112  +  /2221 
+  167r4{(2/3)r3(  -14-  2c3)}2M?/00o2  4-  AII0022  4-  ^3/0222 
4-  A  4(  /0020  4-  2/1120  —  2/1120)  —  2/I1/I2/0012  +  2/1 1/4  3  /0012 
4-  2/4i /44( /ooi  1  /mi)  2/42/43  4)122  2/42/44( /0021 

—  /1 1 2 1 )  4"  ^3/44( /0121  —  /1 22 1 ) I  ~  ((3/ 16)L,L>,/37r4  j2[/4 5/0002 

4-  A 6/0022  4"  "4  7/0222  4~  /4|(/oo2o  2/u2o  4-  /2220) 

4"  2/45/46/0oi2  4“  2/4 5/4 7 An  12  4"  2/4s/48(A)oii  4"  Am) 

4-  2/4^ /47 4)i 22  4-  2/46/48(/oo2i  —  A121)  4“  2/47/48(/oi2i  —  A221)] 

+  { ( 10/3)f 3  }7r2[/45(4)oi2  -  A  112)  "46(4)022  ~  A 122) 

4-  /47( /01 22  A222)  4-  /48( /0021  2/u2i  4-  A221  )](/Vx/A?  -  NJL2)2 

—  (/0022  2/1122  4"  A222X A^/A2  —  Ny/ L2  =  0  (A-3-1) 

where 

/4,  =  sin9  cos9(sin2  0/A?  +  cos2  0/A2) 

/42  =  sin9  cos9(cos2  0/A;J  4-  sin2  0/A2)£u 

A3  =  sin9  cos9(cos2  0  —  sin20)(l/A?  —  1/L2)£n 

/44  =  2  sin9  cos9(cos2  0  —  sin2  0)(1/A2  —  1/A2) 

/45  =  sin4  0/A?  +  6au>26  cos2  0/AxA2  +  cos4  0/A2 

A6  =  (sin4  0/A4  4-  6oiv2d  cos2  6/L2L2  4-  cos4  0/A2)£n 

/47  =  2£,,  cos2  0  sin2  0/A?  -I-  2(cos4  0  4-  sin4  0  —  4  cos2  0  sin2  0) 

X  EJ(LxLy)2  4-2 cos2  0  sin2  0/A4 


Elevated-Temperature  Creep  Buckling  of  Thermoplastic-Matrix  Composites  207 

Aa  =  4  cos2  $  sin2  6/ Li  +  4(cos4  $  4  sin4  6)/(LxLy)2  +  4  cos2  $  sin2  $/ £4 

(A-3-2) 

The  expanded  form  for  Equation  (2-3-2b)  is  shown  as  Equation  (2-3-4).  The 
explicit  form  of  5,  to  BIS  are: 

B\  =  sin8  $ 1 0002  4  4  sin4  $  cos4  ${£11/0222  4"  4£n(/oi2i  4  /mi)  4  4(/0o2o 

—  2/1120  4  /2220)  1  4  £n  cos8  d I Q022  4  4  cos2  $  sin6  $[£11/0112  4  2/oon 

-  2/,,,,]  4  4 £, ,  cos6  $  sin2  $[£,,/om  4  2  /002i  -  2 /, ,2i  J 
4  2EX  cos4  sin4  $/00 12 

B2  =  sin4  $(sin4  $  4  cos4  $) £t  1/0112  4  sin6  $  cos2  $[£i,/ooi2  4  /00 02 

—  4/(/oou  —  /nu)]  4  2  cos2  $  sin2  $(sin4  $  4  cos4  $)£,,[ £11/0022 

4  2 ( Iq  121  /122 1 )]  4  2  sin4  $  cos4  $[£11/0122  4  £11/0112  4£n(/0i2i 

—  /l22l)  4  2£n(/o021  ~  /ll2l)  4  2(/oo„  —  /llll)  “  8(/0020  —  2/h20 

4  £220)]  4  cos4  $(sin4  $  4  cos4  $(£2,/0, 22  4  sin2  $  cos6  $£,,[£,1/0o22 
4  /0012  4(/0021  /l  1 2 1 )  ] 

B3  =  cos2  $  sin6  $(£,,/ooi2  4  /00 02  -  2 IE  ,/0,,2)  4  sin4  $(cos2  0 
sin2  $)2( /oon  4  /nn)  4  2  cos4  $  sin4  $[£11/0122  4  £11/0112 
2  £  n  /0222  4  2£n(/oo2i  —  /1 1 2 1 )  4  2(/oon  /nn)  4£ti(/oi2i 

-  /1221)]  4  2  cos2  $  sin2  $(cos2  $  -  sin2  $)2[ £n( /0121  -  A221) 

4  2( I 0020  ~  2 /]  1 20  4  £220)]  4  cos6  $  sin2  $£n(£n/o 022  4  /0012 

-  2£,,/0,22)  4  £,,  cos4  $(cos2  $  -  sin2  $)2( /002.  -  Am) 

=  sin8  $ £,, /0012  4  2  cos2  $  sin6  $( £, , /0012  4  2/00n  -  2/,,,,) 

4  sin4  $  cos4  $ l 000 2  4  2  cos2  $  sin6  $£,,(£,1/0.22  4  2/002.  -  2/, ,2,) 

4  4  sin4  $  cos4  $[£2/o222  4  4£x(/0i2,  —  /mi)  +  4(/oo2o  —  2/n2o 
4  / 2220)  1  4  2  cos6  $  sin2  $(£,,/on2  4  2/00n  -  2/,,,,) 
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4-  sin4  0  cos4  dEuIoo22  4-  2  cos6  0  sin2  0£11(E,,/0, 22  +  2 /002 1  -  2/im) 
"4“  COS  0£n£oi2 

=  (sin4  0  +  cos4  d)E2uI0222  4-  2  sin2  0  cos2  0(sin4  0  4-  cos4  0)1  £(1/0122 
+  £11/0112  4£u(/0i2i  —  £221)  4-  sin4  0  cos4  0[ £(, /0022  4*  /0002 

16(/0o2o  —  2/  1120  4"  £220)2  £1 1 /0012  —  8(  /001 1  /ill!  )  8  £j  1  ( / 002 1 

£121)] 

B6  =  4-2£,,  sm2  0  cos2  0(sin4  0  4-  cos4  0)(£,,/o  122  4-  /0, 12  -  2£,,/o222) 

4-  £u(sin2  0  —  cos2  0)2(sin4  0  4-  cos4  0)( /01 21  —  £221) 

4"  sin  0  cos  0[ £  1 1  /0022  4~  2 £n £012  4“  /0002  2£(i£i22  —  2£uoi22 

4£ll(£o2l  £l2l)  4(£oil  /1111)  4"  8£,,(/0,2,  / 1  22  1 )  J 

4-  cos2  0  sin2  0(cos2  0  -  sin2  0)2[£,,(/Oo2i  -  £.2.)  4-  ( /00, ,  -  £1.,) 

~  4(4,020  ~  2/u20  4*  £220)] 

B7  =  4-sin4  0  cos4  0[£(,/oo22  4-  2£ll/00,2  4-  £002  4-  4£f , /0222 

—  4 £1 1/0112  ~  4£2 1 /0, 22 ]  4-  cos2  0  sin2  0(cos2  0  -  sin2  0)2[£,,( /002, 

/1121)  4"  (£1011  /mi)  2 £1 1  ( £  121  £221)]  4-  (cos2  0 

sin2  0)4(/oo2o  2  £  1 20  4-  £220) 

=  sin4  0(sin4  0  4-  cos4  0)£(,£, 22  +  £.,  sin6  0  cos2  0[£,,/oo22  4-  /00!2 

—  4(/002i  -  £m)]  4-  2  cos2  0  sin2  0(sin4  0  -I-  cos4  0) £,,[£, , /0222 
4-  2 ( /0 , 2 1  -  £221)]  4-  2  sin4  0  cos4  0[£(,/o,22  4-  £,,/0,, 2 

—  4£„(/o,2i  —  £221)  4-  2£,,(/oo2.  —  /1121)  4-  2(/oon  —  /nil)—  8(/0020 

—  2/1120  4-  £220)]  +  cos4  0(sin4  0  4-  cos4  0)£u £112 
4"  sin  0  COS  0 [  £1 1  /00 1 2  4“  £002  4(/ooii  —  /nil)] 

^9  cos2  0  sin6  0£i i( £1 1  £002  4"  £012  —  2/£,,£,22)  4"  sin4  0(cos20 
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-  sin2  9)2En(IQ02l  -  I i\ii)  4-  2  cos4  0  sin4  9[E\J0l22  +  £11/0012 
2£i, 4)222  4"  2  £,  i  (  /o02 1  ~  A  1 2 1 )  +  2(/ooi!  —  /nil)  —  4£,,(/0, 21 

-  A 22 1 ) ]  +  2  cos2  6  sin2  0(cos2  0  -  sin2  0)2[£,,(/o,2,  -  /mi) 

4"  2(  /0020  2/u2o  4~  /2220)]  cos6  0  sin2  0( £11/0012  4"  /0002 

-  2 E i , /01 1 2)  4-  cos4  0(cos2  0  -  sin2  0)2(/oon  -  Am) 

510  =  sin8  0  £2 , 1 0022  4-  4  sin4  0  cos4  0[£A/O222  4-  4£,,(/0,2i  -  A221) 

4"  4(  /002o  —  2 /,  120  4-  A220  )J  4-  cos8  0/0022  4-  4  cos2  0  sin6  0(  £f ,  /0i 22 
4-  2£,,(/0021  /i  1  2  1  )]  4-  4  cos6  0  sin2  0[£r/o,,2  +  2/oon  -  2/,n] 

+  2£n  cos4  0  sin4  0/ool2 

5,,  =  (l/£4)[sin4  0(/Ooi2  -  Am)  4-  2  cos2  0  sin2  0£11(/O1,2  -  A. 22) 

4-  4  cos2  0  sin2  0(/OO2.  -  2/,  121  4-  A221)  4-  £„  cos4  0(/oo 22  -  A. 22)] 

4-  [8/(£J£2)][(sin4  0  +  cos4  0)£u(/o,22  -  A222) 

2  cos2  0  sin2  0 £, i ( /oi  1 2  —  /1122)  —  4  cos2  0  sin2  0(/oo2i  —  2/,, 2 , 

+  A221  )]  +  ll6/(A2£2)][cos2  0  sin2  0£n(/Oo22  -  / 1122) 

+  cos2  0  sin2  0(/o,,2  -  A. 22)  ~  2  cos2  0  sin2  0 £, ,( /0122  ~  A222) 

4-  (cos2  0  -  sin2  0)2(/Oo2.  -  2 / , , 2 ;  +  A22.)]  +  ( 16/L4)[sin4  0£11(/OO22 

-  A122)  4-  cos4  0(/Ooi2  -  A112)  +  2  cos2  0  sin2  0£,,(/o,22  —  A222) 

4-  4  cos2  0  sin2  0(/OO2i  -  2 /,  12,  4-  A221)] 

5,2  =  (16/£4)[sin4  0(/oo  12  -  A112)  +  2  cos2  0  sin2  0£, ,( /0, 12  -  Am) 

4-  4  cos2  0  sin2  0(/OO2>  -  2/, ,2I  +  A22O  4-  £,,  cos4  0(/OO22  -  A, 22)! 

4-  [8/(£2£2)][(sin4  0  +  cos4  0)£,,( /OI22  -  /1222) 

4-  2  cos2  0  sin2  0£i,(/om  Am)  —  4  cos2  0  sin2  0(/oo2i  —  2/,,2i 
4-  A221)]  +  [16/(A2A2))[cos2  0  sin2  0£,,(/Oo22  -  A122) 
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4-  cos2  0  sin2  0(/ou2  -  /1122)  -  2  cos2  0  sin2  0Ai,(4m2  -  /1222) 

4-  (cos2  0  -  sin2  0)2(/OO2i  -  2/u2i  +  /222O]  +  (1/A4)[sin4  0A,,(/OO22 

—  Am)  4-  cos4  0(/Ooi2  -  Am)  +  2  cos2  0  sin2  0A,,(/O, 22  -  A222) 

4-  4  cos2  0  sin2  0( I 002 1  -  2/,i2i  +  A221)] 

Bl3  =  7o022  ~  2/il22  +  A222 

A,4  =  cos2  0  sin6  0  /0002  +  cos6  0  sin2  0AA/oo 22  +  cos2  0  sin2  0(cos2  0 

sin2  0)2(Af, /0222  "H  4£,i(/o,2,  A221)  4-  4(/0o2o  ~  2/n2 0  4"  A22o)J 

—  2  cos4  0  sin4  0£i  1/0012  —  2  cos4  0  sin2  0(cis2  0  —  sin2  0)A,,[  Ai,/0i22 
4-  2(/002i  -  Am)]  +  2  cos2  0  sin4  0(cis2  0  -  sin2  0)A,,[ An/om 

+  2(/0021  —  A 121)] 

A, s  =  cos6  0  sin2  0  /0002  +  cos2  0  sin6  0A?,/OO22  +  cos2  0  sin2  0(cos2  0 

—  sin2  0)2[ A2, /0222  +  4An(/oi2i  —  A221)  +  4(/002o  ~  2A,2o  +  A220)] 

—  2  cos4  0  sin4  0A,,/Ooi2  —  2  cos2  0  sin4  0(cis2  0  —  sin2  0)A,,[ Au/0,22 
+  2(/002i  —  Am)]  +  2  cos4  0  sin2  0(cis2  0  —  sin2  0) A, , ( Eul0i22 

+  2(  I 002 1  A121)] 

A, 6  =  cos4  0  sin4  0/0002  —  cos2  0  sin2  0(cos2  0  4-  sin2  0)/oo  12 

4-  cos4  0  sin4  0  A,  1/0022  +  (cos3  0  sin  0  —  cos  0  sin3  0)cos  0  sin  0(cos2  0 
sin2  0)[ A2 1  /0122  4-  2An(/002i  —  A121)  4“  Au/0112  4-  2(/oon  /nil)] 

—  cos2  0  sin2  0(cos2  0  -  sin2  0)2[ Af , /0222  +  4All(/0l2i  -  /l221) 

4"  4(4)020  2  /,  1 20  4-  A220)] 

A, 7  =  (4/Aj)[cos  0  sin3  0(/OOi2  -  A 112)  -  cos3  0  sin  BEU(10021  -  I)i22) 

4-  cos  0  sin  0(cos2  0  -  sin2  d)Eu(I0i21  -  A222)  4-  2  cos  0  sin  0(cos2  0 
sin2  0)( /002 '  27,12,  4"  / 2221 ) ]( 1/A2)[sin  0  cos3  0(/ooi2  —  /1112) 
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-  sin3  0  cos  0£\i(/Oo22  —  Aim)  —  cos  0  sin  0(cos2 j  -  sin2  j)Ex flom 

-  )  -  2  cos  0  sin  0(cos2  0  -  sin2  0)(/OO2 1  -  2Iu2l  +  I2 22,)] 

Bls  =  (1/ Z/2)[cos  0  sin3  0(/OOi2  -  Am)  -  cos3  0  sin  0£,,(/Oo 22  -  A 122) 

4-  cos  0  sin  0(cos2  0  —  sin2  6)En(I0l22  —  A222)  +  2  cos  0  sin  0(cos2  0 

-  sin2  0)(/oo2i  -  2 A  12i  +  /222i)](4/L2)[sin  0  cos3  0(/oo.2  -  A, 12) 

-  sin3  0  cos  0£ii(/oo22  —  A122)  —  cos  0  sin  0(cos2  j  —  sin2  j  )£n(/0i22 

-  A222)  -  2  cos  0  sin  0(cos2  0  -  sin2  0)(/oo 21  -  2  A 121  4-  A221)] 

(A-3-3) 


APPENDIX  4 

The  expanded  form  of  Equation  (2-3-2a)  for  the  cross-ply  laminate  is: 


47T4[P?/ooo2  4"  P 2/0022  4-  P  3/0222  4"  P  \{Iq020  2/  1120  4*  /2220)  4-  2PiP 2/0012 


4-  2P1P3/0112  4-  2/>1/>4(/oou  —  Am)  4"  2P2P3/0122  4-  2P2P4(I0o2x 
-  A121)  4-  2P3P4(Iox2x  -  A22.)]  -  5tt2[(/Vx/ZJ)  -  (/V,/L2)][P,(/0o,2 


/1U2)  4"  P'l(  Zo022  /l  1 22)  4"  P3  (  A 1  2  2  A222)  4"  P  /q02  1  2/1121 


4”  A221)]  4-  [(/Vx/£2)  —  (AA/L2)]2(/0o22  —  2/1122  4-  A222) 

[64 Axy/(97r 2LxLy )] 2( /0022  2/1122  4-  A222)  =  0  (A-4-1) 

where 


P,  =  2r3/(3Lx4)  4-  14r3/(3L4) 

P2  =  14r3/(3L4)  4-  2t3/(3Ly) 

P3  =  32Elit3/(3L2xL2) 

P4  =  64t3/(3L2L2)  (A-4-2) 

The  expanded  form  of  Equation  (2-3-2b)  for  the  cross  ply  laminate  is: 

^  [0l05/oOO2  4"  Q4Q&I0022  4-  Q2QbIo222  4"  A)020  —  2  /,  110  +  ^22  20  ) 

+  <e,e8  +  c,es)/oo,2  +  (e,e6  +  qsq,)i<,„,  +  «2,e. 
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+  QiQi){IoOl\  —  h  111)  +  (0406  +  0802)^0122  +  (0407 
4"  0803)(^OO21  —  h  1 2 1 )  4*  (0207  +  060a)(  A)  122  —  A222)] 

-  ir2(Wx/L2)[(40,  +  0s)(/oo,2  -  /.III)  +  (404  4-  08)(/oO22 

—  /1122)  +  (402  +  06 )  ( -A)  1 2  2  ~  /1222)  +  (403  +  07)(/oO21 
~  2/ii2,  4"  /2221)]  —  T2(N,/L2r)l(Q:  +  4es)(/<K,,i  -  /lll2) 

+  (04  +  4^8)(/o022  —  A122)  +  (Q,  +  4G«)(/o.22  -  /l222)  +  (03 
+  40,)(/oc2,  -  2/, ,2,  +  /222,)1  +  [(4Af,!/LJ  +  nN,N,l(LlL}) 

+  4 NllLiKIoou  ~  2 /1122  +  /222)  +  (128W„(97riL,2.,)]J(/0022 

2/1122  +  I2222)  =  0  (A-4-3) 

where 

=  (2t3/3)(l/Ll  +  1  \2/Ly) 

02  =  128f3£„/(3L2L2) 

(3a  =  256/3/(3£x2£2) 

04  =  (2t3/3)Eu(l/Lt  +  16/L4) 

(2s  =  ( 2t3/3)(\6/Ll  +  7/£4) 

06  =  \2St3  ElJ(3L2xLy) 

07  =  256/3/(3L2£2) 

08  =  (2r3/3)£n(7/£?  +  16/£y4)  (A-4-4) 

APPENDIX  5 


ir2£,£y/  2jr  A,  /Vy\ 

-4~(T  -  zj  -  zj)  - 

32  / 

-  —  (2ir25M  4-  Nxy)  ?r2£x£y 


(A-5-1) 
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-  y  (2t r25u 


A, 

L2 


+  A^) 


'T?)~T  (2v‘s“  +  Nx,) 

*‘L,Lyl  4N,  Af,\ 

TT-Ij) 


Cl2 


c2, 


=  0  (A-5-2) 


Using  either  one  of  two  equations  in  above  each  equation,  the  mode  shape  ratio 
can  be  obtained. 

Equation  to  solve  for  the  time-dependent  mode  ratio  cl2/c2i  for  the  angle-ply 
laminates  is: 

[97r2L*Ly/128][(7r2/L;J){sin4  0/oo ,2  4-  2  cos2  0  sin2  6EXXI0 x22  +  4  cos2  0  sin2  0 
X  (/002 1  -  Am)  +  Eu  cos4  flW)  4-  {87r2/(L2L2)j{£u(sin4  0 
4-  cos4  0)/Oi22  4-  sin2  0  cos  0(£i  1  /0022  4“  /00,2  ~  4/oo2,  4-  4/21 121)} 

4-  {167r2/(L2L2)j{sin2  0  cos2  d(EuIo022  4-  /0012  —  2EXXI0X22) 

4-  (cos2  0  -  sin2  0)2(/Oo2,  -  /i,2,)ll  167r2/(L4)j{sin4  6ExxIOQ22 
4-  2  cos2  0  sin2  6ExxI0i22  4-  4  cos2  0  sin2  0(/OO2,  —  /1121) 

+  cos4  0/ooi 2 1  -  {3/(2 t3))(NJLi  +  4 Ny/Lt)(Ioo22  -  /„22)]c12 
-I-  (2/3)7r2[(l/L^){cos  0  sin3  0 /00 1 2  —  cos3  0  sin  6EXXI0 022 
4-  cos  0  sin  0(cos2  0  —  sin2  d)(ExxIox22  4-  2/002, 

—  2/u2i)}(1/L2){cos  0  sin3  0/oo,2  —  cos3  0  sin  0£,,/o 022 

—  cos  0  sin  0(cos2  0  -  sin2  0)(En/Oi22  4-  2/002,  -  2/,,2,)l]c2, 

—  [3/(2r3)]/V,y(/0022  ~  / 1 1 22)^2 1  =  0  (A-5-1) 

Equation  for  the  time-dependent  mode  ratio  cxxlc22  of  the  cross-ply  lami¬ 
nates  is: 

‘K2LxLyl4\p2\P ,(/ooi2  —  /,  1 1 2 1  +  P 2U 0022  ~  /, , 22)  +  Pi(I\ 0,22  —  /1222) 

+  P*(Io02l  ~  2/1121  "b  /2221)}  (A /,/E2  4"  A/y/L2)(/oo22 

2/1122  4"  /2222)]Cii  (32/9)/VIy(/0022  2  /t  1 22 

4-  /2222)c22  =:  0 


(A-5-2) 
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Equation  for  the  mode  ratio  cx2/c2x  of  the  cross-ply  laminates  is: 

TC2LxLy/A'[p2[(4Qx  +  0s)(/()O12  ~  Ain)  +  (404  +  0s)(/oo22  —  A122)  +  (402 
+  08)(/oi22  ~  A222)  +  (403  +  07)(/oO21  —  2  \X2X  +  /222l)l  —  (NXIL\ 

+  4Ny/Ly2)(/0022  —  2/i|22  +  /2222)]Cn  ~  (32/9)  N xy  (  /q022  —  2Al22 
+  l22i2)c2x  =  0  (A-5-3) 

where  the  constants  P,  and  0,  are  the  same  as  those  in  Appendix  4. 

5.  REFERENCES 

1.  Hashin,  Z.  “Analysis  of  Properties  of  Fiber  Composites  with  Anisotropic  Constituents,”  Journal 
of  Applied  Mechanics,  46:543  (1979). 

2.  Massonet,  Ch.  “Buckling  Behavior  of  Imperfect  Elastic  and  Linearly  Viscoelastic  Structures" 
International  Journal  of  Solids  and  Structures,  10:755  (1974). 

3.  Lekhnitskii,  S.  G.  Anisotropic  Plates.  New  York,  NY.Gorden  and  Breach  Science  Publishers, 
pp.  382,  437  (1968). 

4.  Williams,  J.  G.  Stress  Analysis  of  Polymers.  London,  UK: Longman  Group  Limited,  p.  148 
(1973). 

5.  Powell,  P.  C.  “Principle  for  Using  Design  Data,”  in  Thermoplastics.  R.  M.  Ogorkiewicz,  ed. 
Cleveland,  OH.CRC  Press,  p.  211  (1969). 

6.  Flugge,  W.  Viscoelasticity.  Waltham,  MA:BlaisdeII  Pub.  Co.,  p.  92  (1967). 

7.  Hilton,  H.  H.  “Creep  Collapse  of  Viscoelastic  Columns  with  Initial  Curvatures,"  Journal  of 
Aeronautical  Sciences,  19:844  (1952). 

8.  Freudenthal,  A.  M.  The  Inelastic  Behavior  of  Engineering  Materials  and  Structures.  New  York. 
NY:John  Wiley  &  Sons,  Inc.,  p.  506  (1950). 

9.  Fiacco,  A.  V.  andG.  P.  McCormick.  Nonlinear  Programming — Sequential  Unconstrained  Mini¬ 
mization  Techniques.  New  York,  NY:John  Wiley  &  Sons  (1968). 

10.  Hoff,  N.  J.  “Buckling  and  Stability,”  Journal  of  the  Royal  Aeronautical  Society,  58(517):  1  (1954). 

11.  Yamamoto,  Y.  “Considerations  of  a  Stability  Criterion  for  Creep  Buckling,”  Journal  of  Mechan¬ 
ics  and  Physics  of  Solids,  18:165  (1973). 

12.  Ross,  A.  D.  “The  Effects  of  Creep  of  Instability  and  Indeterminacy  Investigated  by  Plastic 
Models,”  The  Journal  of  the  Institution  of  Structural  Engineers,  24:421  (August  1946). 

13.  Vinogradov,  A.  M.  “Buckling  of  Viscoelastic  Beam-Columns,"  27th  Structural  Dynamics  and 
Materials  Conference,  American  Institute  of  Aeronautics  and  Astronautics  Inc.,  1:502  (1986). 

14.  Booker,  J.  R.  et  al.  “Stability  of  Viscoelastic  Structural  Members,"  Civil  Engineering 
Transactions,  p.  45  (1974). 

15.  Booker,  J.  R.  “A  Method  of  Solution  for  the  Creep  Buckling  of  Structural  Members  of  a  Linear 
Viscoelastic  Material,"  Journal  of  Engineering  Mathematics,  7(2):  101  (April  1973). 

16.  Lin,  T.  H.  “Stress  in  Columns  with  Time  Dependent  Elasticity,”  Proceedings  of  the  1st  Midwest¬ 
ern  Conference  on  Solid  Mechanics,  p.  196  (1953). 

17.  Miyase,  A.,  A.  W.-L.  Chen.  P.  H.  Geil  and  S.  S.  Wang.  “Anelastic  Deformation  and  Fracture  of 
Thermoplastic-Matrix  Fiber  Composite  at  Elevated  Temperature,  Part  II  Undirectional  Compos¬ 
ite  Laminate,”  American  Chemical  Society  20th  Great  Lake  Regional  Meeting,  Milwaukee  (June 
1986). 

18.  Froberg,  C.-E.  Introduction  to  Numerical  Analysis.  2nd  ed.  Reading.  MA.  Addison-Wesley  Pub¬ 
lishing  Company,  p.  348  (1969). 


Elevated-Temperature  Creep  Buckling  of  Thermoplastic-Matrix  Composites 


215 


19.  Ferry,  J.  D,  Viscoelastic  Properties  of  Polymers.  New  York,  NY,  section  11  (1970). 

20.  Wang,  S.  S.,  E.  C.  Klang,  A.  Miyase  and  A.  Dasgupta.  “Viscoelastic  Shear  Behavior  of  the 
Thermoplastic  Matrix  Composite  at  Elevated  Temperatures,”  ASME,  Winter  Annual  Meeting, 
Boston  (December  1987). 

21.  Wilson,  D.  W.  and  J.  R.  Vinson.  “Viscoelastic  Buckling  Analysis  of  Laminated  Composite  Col¬ 
umns,”  Recent  Advances  in  Composites  in  the  United  States  and  Japan ,  ASTM  STP  864,  J.  R. 
Vinson  and  M.  Taya,  editors,  American  Society  for  Testing  and  Materials,  Philadelphia,  p,  368 
(1985). 

22.  Wang,  S.  S.,  A.  Miyase,  P.  H.  Geil  and  A.  Chen.  “Anelastic  Deformation  and  Fracture  of  Ther¬ 
moplastic  Matrix  Composites  at  Elevated  Temperatures”  (1988). 


